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Phonons in single crystals of PrFeAsOi_ H are investigated using high-resolution inelastic x-ray 
scattering and ab-initio pseudo-potential calculations. Extensive measurements of several samples 
(y ~0, 0.1 and 0.3) at temperatures spanning the magnetic ordering temperature (TV ~145 K for 
y ~0) and the superconducting transition temperature (T c — 36 K for y ~0.1 and T c — 45 K for 
y ~0.3) show that there are some changes in phonon spectra with temperature and/or doping. We 
compare our measurements with several ab initio pseudo-potential models (nonmagnetic tetragonal, 
oxygen-deficient 7 / g supercell, magnetic orthorhombic, and magnetic tetragonal) and find that the 
experimentally observed changes are much smaller than the differences between the experimental 
data and the calculations. Agreement is improved if magnetism is included in the calculations via 
the local spin density approximation, as the Fe atomic motions parallel to the ferromagnetic order- 
ing direction are softened. However, the antiferromagnetically polarized modes remain hard, and in 
disagreement with the experimental data. In fact, given the increasing evidence for anisotropy in the 
iron pnictide materials, the phonon response is surprisingly isotropic. We consider several modifica- 
tions of the ab initio calculations to improve the agreement with the experimental data. Improved 
agreement is found by setting the matrix to zero (clipping the bond) between nearest-neighbor an- 
tiferromagnetically aligned Fe atoms in the magnetic calculation, or by softening only the in-plane 
nearest-neighbor Fe-As force constant in the nonmagnetic calculation. We discuss these results in 
the context of other measurements, especially of phonons, for several FeAs systems. Fluctuating 
magnetism may be a partial explanation for the failure of the calculations, but seems incomplete 
in the face of the similarity of the measured phonon response in all the systems investigated here 
including those known to have static magnetism. 



PACS numbers: 74.25. Kc, 74.70.Xa, 78.70.Ck 
I. INTRODUCTION 



Since the discovery of superconductivity in 
LaFeAs(0,F) at 26 K (Ref. [l|), iron-arsenides and 
related compounds have been the subject of enormous 
scientific attention. These compounds have a layered 
structure consisting of iron and pnictide/chalcogenide 
atoms ( "Fe-As layer" ) with each iron atom surrounded 
by a tetrahedron of pnictide/chalcogenide atoms (c.f. 
Fig. [1]). The iron-based superconductors have been 
classified into, mainly, four groups according to the 
kind of buffer between the Fe-As layers: i?FeAsA (R: 
rare earth, X = O or F) with ZrCuSiAs type structure 
("1111" ^-I, i?Fe 2 As 2 (R = Ba, Sr, Ca) with ThCr 2 Si 2 
type structure ("122" )S&, ^4FeAs (A = Li, Na) with 
Cu 2 Sb type tetragonal structure ("111") 10 ! 11 , and 
a-FeCh (Ch = S, Se, Te) with a-PbO type structure 
("ll") 12 ' 13 . Doping electron or hole carriers into the 
Fe-As layers, or application of pressure (external or 
chemical) results in superconductivity at low tem- 
perature. Some of the "1111" type materials show 
superconductivity above 50 K (Refs. 0-0), which is the 



highest observed superconducting transition tempera- 
ture (T c ) outside of the copper oxide family of high-T c 
superconductors. Recently more complex compounds 
having a buffer layer of perovskite structure have been 
synthesized pursuing a higher T c (RefsJ^r— ). However, 
the highest T c is still achieved in a "1111" material. 

It has been suggested by many groups that the phys- 
ical properties of iron-arsenic compounds are very sen- 
sitive to the crystal structure, especially to the pnic- 
tide/chalcogenide height above the Fe plane and/or 
the angle of the As-Fe-As bonds 2 ^— . This is fur- 
ther emphasized by the possibility to dramatically af- 
fect the superconducting properties of some material by 
applying pressure, with T c increased in LaFcAs(0,F) 
from 26 to 43 K under an applied pressure of 4 GPa 
(Ref. [H), and in Fci.oiSe from 8.5 to 36.7 K under 
a pressure of 8.9 GPa (Ref. l27t ). This sensitivity sug- 
gests the possibility of large coupling between lattice 
motion and electronic or, perhaps, magnetic structure. 
Meanwhile, early first principles calculations of phonon 
properties suggested that the superconductivity in iron- 
arsenic compounds is not phonon mediated 2 ^—. How- 
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TABLE I: PrFeAsOi-y samples. The thickness of the samples 
was a few tens pm. The sample 'doped- 1' was also used in a 
previous study fRef. l37h. 



FIG. 1: (Color online) Crystal structure of _ReFeAsO {Re: 
rare earth). The arrows on Fe atoms show the magnetically 
ordered pattern observed for T < Tn (Refs. QHJ). Dashed 
(red) arrows indicate the direction of ferromagnetic (F) [an- 
tiferromagnetic (AF)] alignment of iron spins in the plane. 



ever, later calculation a 31 ' 32 demonstrated the strong de- 
pendence of the iron magnetic moment on the pnic- 
tide/chalcogenide atomic position, suggesting some in- 
terplay between phonons and magnetism, while others 
have noted increased electron-phonon in magnetic calcu- 
lations by factors of -2 (Ref. to -10 (Ref. 

Experimentally phonons have been measured using 
various methods. Inelastic neutron and/or x-ray scatter- 
ing (INS/DCS) is a powerful tool to investigate phonon 
structures, and many measurements have been already 
reported 3 -^—. Many of these reported small changes with 
doping and/or temperature. The experimental phonon 
data are compared to first principles calculations, and 
the introduction of magnetic effect is thought to be 
important^—. We will discuss them in detail later. 
Here we note that with the exception of Refs. and 
|49| all single crystal work has focused on the "122" sam- 
ples, despite their lower T c , because the "1111" materials 
remain relatively difficult to grow. 

In the present paper, we report phonon dispersion 
measurements along various symmetry directions mea- 
sured on single crystals of PrFeAsOi_ a using IXS. This 
is a "1111" compound that shows relatively high T c up 
to 45 K for the samples discussed here and 49 K for 
the Pr family. We also perform several first principles 
calculations, and compare them carefully with the ex- 
perimental data. Our detailed observation confirmed 
the isotropy of phonon dispersion even in the antifcrro- 
magnetically ordered phase compared to the first princi- 
ples calculation considering the magnetic effect. Though 
some similar observations, calculations, and the compar- 
isons have been already discussed for other iron-arsenide 
compounds 3 ^—, the similarity of our results to other 
works ensures that the present work contains the gen- 
eral features of iron-arsenides. 

We then consider several modifications of the first prin- 
ciples calculations to better interpret our results. Two 
modified models are found to have good agreement with 
the experimental data. Interestingly, one is the model 
that adds additional in-plane anisotropy, while the other 
preserves the isotropy of our nonmagnetic first principles 
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calculation. 

This paper is organized as follows. The details of sam- 
ples investigated, IXS measurements, and first principles 
calculations are provided in Sec. [TXJ The experimental 
results and the comparison with calculations are given 
in Sec. IIIII Several modification model of first principles 
calculations are detailed in Sec. [TV] Section [V] is dedi- 
cated to the discussion of the results and a summary is 
given in Sec. EH 



II. EXPERIMENTS AND CALCULATIONS 
A. Samples 

The investigated PrFcAsOi_ y single crystals arc sum- 
marized in Table HI They were prepared using high- 
pressure growth as described in Ref. [52- The first sam- 
ples were relatively small — 0.1 x 0.15 x 0.02 mm 3 , but 
with improvements in growth techniques, they became 
comfortably large, —0.5 mm in the ab plane. Sample 
thickness varied from about 0.02 to 0.05 mm. X-ray 
diffraction on a four-circle diffractometer was used to ver- 
ify that all samples were single grains, with, typically, a 
mosaic spread of about 1°. The T c 's of the studied su- 
perconducting samples were found to be 36 and 45 K by 
measuring magnetic susceptibility using a superconduct- 
ing quantum interference device (SQUID) magnetometer. 
The electrical resistivity in the parent sample showed 
an abrupt change at —149 K and the derivative had a 
maximum at —139 K (Fig[2]). According to the previ- 
ous transport measurements in iron pnictide o 53 ' 54 , these 
temperatures correspond to the structural and magnetic 
transition, respectively. These temperatures are slightly 
higher than the values given in Refs. [H and [5^, possibly 
because we are more nearly at the precise stoichiometric 
composition. There is a small anomaly in p around 12 
K where the Pr spins order. This is nearly the same as 
previously reported 3 -^. 

We use tetragonal notation for all samples, with axes 
along the next-nearest-neighbor iron atoms. In the par- 
ent sample the lattice becomes orthorhombic below T s , 
followed by ordering of the Fe magnetic moments be- 
low Tjv- When required, we take the ferromagnetically 
(antifcrromagnetically) ordered direction in the Fe plane 
to be parallel to <110> (<110>) as shown in Fig. [T] 
The lattice constant along <110> or the ferromagncti- 
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FIG. 2: (Color online) Electric resistivity (p) and the deriva- 
tive (dp/dT) of a PrFeAsO sample near the structural and 
magnetic transition temperatures. Solid (dashed) lines corre- 
spond to the measurements with decreasing (increasing) tem- 
perature. Inset shows a larger temperature range for p. The 
deviation in p from slow increase with decreasing temperature 
or an abrupt increase in dp/dT corresponds to the structural 
transition temperature T s from tetragonal to orthorhombic 
symmetry, while the inflection point in p or the maximum 
in dp/dT corresponds to the magnetic transition temperature 
Tjv. A small upturn below 12 K indicated by an arrow in the 
inset corresponds to the ordering of Pr magnetic moments. 



or <110> directions approximately vertical. 

The use of the two-dimensional analyzer array allowed 
parallelization of measurements for either nearly pure 
longitudinal or transverse modes dispersing out from the 
r points. For example, with <010> vertical and the 
<100> in the horizontal scattering plane, we can mea- 
sure simultaneously four Q points along approximately 
<100> longitudinal direction or three Q points along ap- 
proximately <010> transverse direction from a (H 0) 
Bragg points using, respectively, a horizontal or verti- 
cal line of analyzers. Likewise, <110> longitudinal and 
transverse modes can be investigated easily with the con- 
dition <110> vertical. The analyzers were not always 
perfectly centered along a symmetry direction, but the 
deviation was typically very small [<0.02 reciprocal lat- 
tice units (r.l.u.) in H and K, and <0.15 r.l.u. in L at 
the momentum transfers we measured]. A typical value 
for the momentum resolution (corresponding to the an- 
alyzer acceptance) is 0.06 r.l.u. in H and K, and 0.12 
r.l.u. in L. 

To analyze the data, the obtained spectra were fit to 
the sum of a resolution-limited clastic peak and several 
Lorentzian phonon peaks: 
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cally ordered direction is shorter than that along <110> 
or the antiferromagnctically ordered direction, a<no> < 
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B. IXS measurements and data analysis 



where Hes(E) is the instrumental resolution function 
measured using elastic scattering from plexiglas, A and 
Bi are constants, Ei and bi are the energy and the width 
of the z-th phonon peak, and n(E) is the Bose thermal 
factor. Fitting to the above function was done using a 
nonlinear Levenbcrg-Marquardt least squares algorithm. 



IXS measurements were performed at BL35XU 
(Ref. 15(1 ) of SPring-8. The x-ray beam from the storage 
ring was first monochromatized to ~eV by a double crys- 
tal, liquid nitrogen cooled, Si(l 1 1) monochromator, and 
then down to 0.8 meV using a Si(ll 11 11) backscatter- 
ing monochromator operating at 21.75 kcV. The incident 
beam was focused, by a bent cylindrical mirror, to a spot 
size of about 70 microns in diameter [full width at half 
maximum (FWHM)] at the sample position. The scat- 
tered radiation was analyzed using 12 spherical crystals 
in a 4 (horizontal) x3(vertical) array on a horizontal 10 
m 29 arm. Each analyzer focused the radiation into a 
separate detector located near the sample. 

The samples were mounted on thin glass rods, which 
were set on a small goniometer head in air for room tem- 
perature measurements, or on the cold finger of a He 
closed-cycle refrigerator for low temperature investiga- 
tions. The penetration depth of PrFeAsOi_ a for 21.75 
keV x rays is about 60 /zm, which, with our thin samples, 
enabled us to measure in a transmission (Lauc) geometry 
for many Q points. The surface normal of all samples was 
along the c axis, and they were aligned with the <010> 



C. First principles calculations 

Several first principles calculations of PrFeAsO i- y and 
LaFeAsO were carried out using VASP— ~— and a pro- 
jector augmented-wave (PAW) metho d 61 ' 62 , in the gen- 
eralized gradient approximation (GGA), and local spin 
density approximation (LSDA). (Table HI)) Phonons were 
calculated via the PHONON package using a direct 
method 6 ^. The supercell size used for each calculation 
is 2a x 2a x c (32 atoms) for "original" , 2a x 2a x c (31 
atoms) for "0 7/8 ", and 2y/2a x 2y/2a x 2c (128 atoms) 
for both "magnetic" and "mag. tetra," where a and c are 
lattice constants of the tetragonal primitive unit cell (cf. 
Table HI)) . Total energies and inter-atomic forces were cal- 
culated for 16, 56, 30, and 30 symmctry-inequivalcnt dis- 
placements for original, O7/8, magnetic and mag. tetra 
calculations, respectively The energy cut-off for plane 
waves was 550 eV, the spacing of k points was less than 
0.1 A -1 and the convergence condition was that the total 
energy difference be less than 1 zxeV. All structures were 
relaxed. The structural parameters are summarized in 
Table [XXJ Stripe antiferromagnetic ordering, as observed 
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TABLE II: Structural parameters from first principles calculations (upper portion) and the experiment (lower portion), "orig- 
inal" is the simplest calculation with nonmagnetic tetragonal structure, while "O7/8" is a supercell calculation with a 12.5% 
ordered oxygen deficiency, "magnetic" and "mag. tetra" are the first principles calculations with antiferromagnetically ordered 
Fe moments (after Ref. @), and the crystal symmetry is orthorhombic and tetragonal (a ort ho = b or tho', Ibam), respectively. 
"O7/8" gives different zp r 's and zas's and the averages are given. In the parent material below Tjv Fe spin moments are aligned 
ferromagnetically along the b or tho direction (||<110>, the notation here as shown in Fig. [T}. 
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experimentally, is stabilized in the magnetic calculations, 
with structural parameters that are nearer to the exper- 
imentally determined values^. 

Magnetic calculations were done using the La pscu- 
dopotential instead of Pr to avoid difficulties in treating 
the spins of the localized Pr / electrons in the core pseu- 
dopotcntial. However, a comparison of the results of non- 
magnetic calculations using La and Pr showed the results 
to be nearly identical for the phonons, with the differ- 
ences between the two calculations being much smaller 
than those between the calculation and the measured 
data. Considering a purely mass difference would sug- 
gest phonon frequency changes of less than 1 percent for 
La/Pr modes. 
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III. RESULTS 
A. IXS measured data 

Figure [3] shows typical IXS spectra of parent PrFeAsO 
and doped PrFeAsOi_ y (doped-2; T c =45 K) at room 
temperature at Q = (3.03 0.06) (near the Brillouin zone 
center, T point) and (3.50 0.00) (Brillouin zone bound- 
ary). An elastic peak and several phonons were observed, 
and fit using the function of Eq. ([TJ . The sum fit curve as 
well as individual phonon lineshapes are shown in Fig. [3] 

Figure @] shows the dispersion relations of PrFeAsOi_ y 
(parent: y =0, doped-1: y —0.1 or doped-2: y —0.3) at 
room temperature along some high-symmetry directions; 
Q = (3+q 0), (3 q 0), (0 q 9), and (q q 9). The area 
of the symbols is proportional to the integrated peak in- 
tensity, Bi, and the error bars give the intrinsic width, 
after subtracting the measured instrumental resolution 
width from |6j/2|. The uncertainty in mode energy is 
<0.5 meV except for some small and/or broad peaks. 



FIG. 3: (Color online) Typical IXS spectra of doped 
PrFeAsOi-y (doped-2; open circle) and parent PrFeAsO 
(closed circle; shifted vertically for clarify) at room tempera- 
ture (a) near T point (Q = (3.03 0.06)) and (b) at Brillouin 
zone boundary (Q = (3.50 0.00)). The lines show the fit to 
the data using Eqn. (QJ. The individual peaks are also shown 
for the doped-2. The arrows below each plot show the peak 
positions based on the fits. Generally the elastic intensity is 
larger near to the T point. 



Though the data in doped PrFcAsOi_ y having slightly 
low T c (doped-1) are plotted only in Fig.@|b), the doped- 
1 and doped-2 samples show exactly the same result at 
Q = (3+q 0) (not shown). 

There is a doping dependence, which is also branch 
dependent. For example, there are three strong branches 
observed at Q = {3+q 0) in Fig. [3(a). With doping, the 
branch dispersing from —32 meV at Q — (3 0) to —27 
meV at (3.5 0) hardens, while the branch dispersing 
from -16 meV at (3 0) to -23 meV at (3.5 0) softens 
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FIG. 4: (Color online) Dispersion relations along some high- 
symmetry directions of doped PrFeAsOi-y (open circle) and 
parent PrFeAsO (closed circle) at room temperature. The Q 
= (3 q 0) data of (b) are from the sample with y ~0.1 (T c — 
36 K; doped-1), and the other ones are from the sample with 
y ~0.3 (T c = 45 K; doped-2). The area of the symbols is 
proportional to the peak intensity [Bi in Eq. {TJ] and the error 
bars are the intrinsic peak widths estimated by subtracting 
the experimental resolution from the fit result. 



slightly. The hardening/ softening may be larger at the 
zone boundary than at the zone center, which can be 
seen more clearly in Fig. [31 However, the differences are 
small (~0.5 meV), and there is no drastic change in the 
overall dispersion between the parent and doped sample 
data. The doping dependence in dispersion relations at 
low temperature below Tjv and T c show features similar 
to Fig. [4] (not shown) . 

Figure [5] shows the temperature dependence of the dis- 
persion in the parent PrFeAsO at Q = (3+q 0), (3 <? 0), 
(0 q 9), and (q q 9). We can see small changes in some 
high energy branches. The mode energy increases with 
decreasing temperature, as may just be the result of ther- 
mal contraction. Measurements about T s (and T/v) near 
Q = (3 0) [Fig. Ufa)], did not show any strong changes 
on crossing the transition temperature, with mode fre- 
quencies being nearly unchanged and most mode intensi- 
ties following the usual Bose factor. That is, the temper- 
ature dependence in doped superconducting PrFeAsOi_ y 
is similar to that in the parent sample, with a slight hard- 
ening as temperature is decreased, and no drastic change 
in the overall dispersion. 

Measurements of the parent PrFeAsO were made along 
the <110> direction at 10 K, well below the measured 
Tn — 139 K, to investigate possible effects of the or- 
thorhombic structural distortion and appearance of mag- 
netism. This direction potentially twins, with the <110> 
and <110> directions becoming distinct. If the beam hits 
a twinned section of the sample, it would be reasonable 
(especially based on the magnetic calculations discussed 
below) to expect phonon splitting as the modes polarized 



FIG. 5: (Color online) Dispersion relations along some sym- 
metry directions for the parent PrFeAsO. The figure style is 
the same as that of Fig. f3] The data at room temperature are 
identical to those in Fig. [4] 
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FIG. 6: (Color online) No magnetic splitting is observed in 
the phonon spectra for PrFeAsO at 10 K. (a) 29 scans near 
(3 3 0)/(3 3 0) Bragg point. The (3 3 0)/(3 3 0) correspond to 
the scattering from antiferromagnetically/ferromagnetically 
ordered (orthorhombic a/6) direction, respectively. (b) 
Phonon spectra at Q — (3.22 3.21 -0.11). The lines are the 
spin polarized ab initio "magnetic" orthorhombic calculation, 
(see Sec. IIII B[) convoluted with the instrumental resolution, 
(c) Dispersion relations. 



in the <110> ferromagnetic ordering direction (shorter 
lattice constant) could have different energies from those 
in the <110> antifcrromagnctic ordering direction. With 
some effort, we were able to find part of the sample where 
twinning was clearly observed, based on the split of the 
(3 3 0) Bragg reflection [Fig. |B{a)], however, even in this 
twinned region, the phonon splitting near ~34 meV ex- 
pected to be ~3 meV from magnetic calculation did not 
appear as shown in Fig. [H^b) . The dispersion is plotted 
in Fig. |Hl[c) for the parent sample from two regions at 
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FIG. 7: (Color online) Comparison of the measured dispersion 
for PrFeAsOi-y (T c = 45 K; doped-2) at room temperature at 
Q = (3+q 0) and various first principles calculations. The 
same data are repeated in each panel for comparison with 
the different calculations: (a) "original" calculation using a 
nonmagnetic tetragonal structure, (b) the "0 7 / g " calculation 
with an oxygen deficiency, (c) the "magnetic" first principles 
calculation, and (d) the tetragonal magnetic calculation. See 
also the text and table [TT| for details. The area of the ex- 
perimental data points (line height of the calculations) show 
the measured (calculated) mode intensity for IXS spectra. For 
clarity, a cutoff has been used in plotting acoustic mode inten- 
sity which becomes large near Bragg points. The low energy 
acoustic mode appearing at (4 — q 0) is probably a transverse 
mode. See text at the end of Sec. IIIIBI 



low temperature as well as those at room temperature. 
The spectral shape, and the dispersion is essentially iden- 
tical. One notes that there is some hardening with de- 
creasing temperature observed in the high energy branch 
[Fig. injc)] and almost no doping dependence in this di- 
rection (not shown). This small doping/temperature ef- 
fects are similar to other phonon modes discussed above 
(Figs. H and |5j}. 



B. Comparison with ab initio calculations 

We now compare the dispersion against the unmodified 
ab initio calculations. Since the effect of carrier doping, 
or sample temperature, on the phonon spectra is rela- 
tively small (see Figs. 0] and O we plot only one data set, 
that from the superconducting PrFcAsOi_j, (dopcd-2). 
Figure [7] shows the comparison between the experimen- 
tal data and the calculations at Q = (3+q 0). Again, 
the area of experimental data points shows the mode in- 
tensity after removal of the Bose factor. For the calcula- 
tions, the line height shows the expected intensity. The 
scale factor (symbol area to line height) is fixed for all 
data sets. The experimental data in Fig. [7] are the same 
as shown previously [doped-2 in Fig. EJa)], but the error 
bars indicating the intrinsic peak width are omitted. 



The agreement between the data and calculation is 
fairly good for low energy branches. However, the branch 
dispersing from ^33 mcV at the T point shows notable 
differences between the different calculations. The mea- 
sured energy is generally significantly lower than that 
of the first principles calculation with a nonmagnetic 
ground state [Fig. [3a) "original"], consistent with our 
earlier work^i. As the superconducting samples are oxy- 
gen deficient, we also made calculations with a deficiency 
(12.5 %) using a large unit cell size 2x2x1 with one 
oxygen removed [Fig.JTJb) "O7/8"]. The agreement with 
the experimental data is similar to the "original" calcu- 
lation, but many added branches appear, as might be 
expected from adding an ordered deficiency. Such addi- 
tional modes are not observed in our measurement, and 
the agreement with the high energy branch in the "orig- 
inal" calculation is not improved. This calculation then 
mostly provides confirmation that the oxygen deficiency 
probably does not have a surprising impact on the de- 
tailed phonon dispersion. This is consistent with the re- 
cent first principles calculation on the doping dependence 
of phonon DOS in LaFeAsOi-^Fj; using the virtual crys- 
tal approximation^. 

We also carried out calculations with antiferromagnet- 
ically ordered Fe moments in an orthorhombic structure, 
essentially that of de la Cruz et al£. The magnetic or- 
der of Pr atoms was not considered as they were substi- 
tuted with La in the calculation. The result is plotted 
in Fig. Etc) (magnetic) . The agreement with the experi- 
mental data is better than the "original" model. We also 
performed a calculation of a tetragonal magnetically or- 
dered material to distinguish between the effects of the 
magnetic moment and crystal symmetry. The result is 
plotted in Fig. [7](d) (mag. tetra), and it is very similar 
to the "magnetic" calculation. This clearly shows the 
phonon softening of the branch at 27-33 meV in the first 
principles calculation comes from allowing the magnetic 
order of the iron atoms: the orthorhombic/tetragonal 
crystal symmetry change has a relatively small effect on 
the calculated phonon dispersion. 

The addition of the magnetism, at first glance, im- 
proves the agreement between the data and the calcula- 
tions. It is worth emphasizing that this is true for both 
the parent material (data not shown) which explicitly 
shows magnetic order, and also the superconducting ma- 
terials [data as shown in Fig.JTJc)] which do not show ev- 
idence of static magnetic order. The latter is somewhat 
surprising, and we will discuss it later again. However, 
while the magnetic calculations do better, they also pre- 
dict splitting of modes that is much larger than that ob- 
served in our data. This is evident in Figs.[7|c) and (d), 
where calculations give a high energy branch (~34 meV) 
that is not observed here. The high energy branch origi- 
nates from the magnetic calculation lifting a degeneracy 
between the ferromagnetically and antiferromagnetically 
polarized modes: the measured phonon energies are con- 
sistent with the energies calculated for modes with Fe 
motions in the ferromagnetic ordering direction [motion 
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FIG. 9: (Color online) Dispersion along Q = (3 + q3 + q0) 
for the parent sample PrFeAsO. In this direction of (c) "mag- 
netic" and (d) "mag. tetra" calculations there are two do- 
mains originated from crystal structure and/or magnetic or- 
der, and the phonon dispersions from each domains are plot- 
ted. See the caption of Fig. [7J for details. The weak, low 
energy, acoustic mode that can be seen in the data is a trans- 
verse acoustic mode (dashed line) that appears in the mea- 
sured spectra, probably due to the finite momentum accep- 
tance of the analyzers. See the text at the end of Sec. IIIIBI 
for discussion. 



FIG. 8: (Color online) The bottom right is a schematic of 
the Fe-As layer: As atoms are located above or below the Fe 
atomic plane as specified. The arrows on Fe atoms specify the 
magnetic spin directions for "magnetic" and "mag. tetra" cal- 
culations, where the ferromagnetic (antiferromagnetic) prop- 
agation vector corresponds to be <110> (<110>) as shown in 
the red dashed arrow F (AF). In (a)-(i) arrows specify atomic 
displacement patterns in an Fe-As layer at several points in 
Figs.[7Jand[I|] Atomic movements along out-of-plane direction 
are specified by "up" or "down" near the atoms. 



A in Fig. [5{a)] but not with the calculated energies of 
modes with Fe motions in the antiferromagnetic direc- 
tion, of which should also appear in the (3 + q 0) direc- 
tion [motion B in Fig. [S^b)]. 

The failure of the magnetic calculations is further con- 
firmed by investigating dispersion along the <110> di- 
rection. A splitting similar to that mentioned above is 
calculated to appear between ferromagnetically polarized 
modes (observed in a longitudinal <110> geometry) and 
antiferromagnetically polarized modes (observed in the 
<110>) in a twinned portion of the parent below TV- 
As discussed previously, Fig. [5] shows the measurement 
of longitudinal modes in the <110> direction, from both 
a single domain and a twinned section of the crystal. The 
phonon spectra are essentially identical, showing no evi- 
dence of the calculated splitting. 

The lack of splitting is also confirmed in Fig. |H] along 
(3 + q 3 + q 0) corresponding to Fig. HJc). In the fig- 
ure, the experimental data are compared with (a) "orig- 



inal" , (b) "0 7/8 " , (c) "magnetic" ,_and (d) "mag. tetra" 
calculations. The <110> and <110> directions are not 
equivalent, as is shown by choosing different colors (gray 
scales). The general features are consistent with those 
of Q = (3+q 0) in Fig. [7J in that the branches at low 
energy are explained fairly well by any calculation, while 
the branches at the energy of 30-35 meV are not. A 
similar feature is also obtained in the data along other 
symmetry directions of Q = (3 q 0), (0 q 9) and (q q 9), 
which arc corresponding to (b)-(d) in Figs. 0] and [5] 

Figure [5] investigates the calculated polarizations re- 
lated to the splitting of ferromagnetically and antiferro- 
magnetically polarized modes in more detail, showing the 
atomic motions at T and M points. For these modes only 
Fe and As atoms move, while rare-earth (La/Pr) and O 
atoms do not. The motions A and B are degenerate in 
the "original" calculation, and split in "magnetic" and 
"mag. tetra" calculations. The pairs D and E, F and 
G, and H and I are the same as the pair A and B. In 
each pair the phonon energy of motions A, D, F, and H 
is lower than that of motions B, E, G, and I. In all cases, 
the motion containing the vibration of Fe atoms along the 
<110> direction (ferromagnetic mode) has lower energy 
than that along the <110> direction (antiferromagnetic 
mode) excepting those modes without in-plane Fe motion 
(i.e., H and I where Fe atoms do not vibrate). 

Finally we note the appearance of a transverse mode 
in the longitudinal spectra, as can be seen in Fig. |H1 
The predominantly transverse character of the mode is 
confirmed by its higher intensity in the "off symmetry" 
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analyzers (data not shown) which have a larger trans- 
verse contribution. Its appearance, in principle, then can 
be explained by the finite momentum acceptance of the 
analyzers [for example AQ ~ 0.029 r.l.u. near the Q = 
(3 0)]. However, the investigation of the calculations 
in both this section, and the next, shows that there is 
also the potential that the mode itself may be of mixed 
polarization. Furthermore, the amount of mixing might 
be some way of selecting between different calculations, 
though such a detailed analysis, relying heavily on both 
the intensity and momentum resolution is beyond the 
scope of the present work. 



IV. COMPARISON WITH MODIFIED MODELS 

The "magnetic" calculations provide the best over-all 
agreement with our data out of those presented above. 
However, there remain significant discrepancies with the 
calculations predicting mode splitting that is not ob- 
served experimentally as discussed above. To better un- 
derstand these continued discrepancies, we consider sev- 
eral direct modifications to the real-space force constant 
matrices resulting from (the interpolation of) the ab ini- 
tio calculations. Since the discrepancies between the cal- 
culations and the data remain generally larger than those 
between different dopings or temperatures, we focus on 
modifying the models to get globally similar character- 
istics to the data. The models are summarized in Ta- 
ble IIIII and Fig. [TUJ For all of them we preserved the 
crystal symmetry, and recalculated the self forces to keep 
the translational invariancc. Optimizations were done by 
hand. 

The simplest model we consider is the "soft Fe-As" 
model of Ref. [3?] where the nearest neighbor Fe-As force 
constant was scaled (reduced) by 30 % [Fig. fTUtaVl] 
from that of the nonmagnetic 'original' calculation. This 
was shown to reproduce the softening of the peak in the 
DOS, but, not completely the dispersion in our previous 
work^i. The comparison with the in-plane longitudinal 
mode along <100> is shown in Fig. [TTTa). While the 
high energy branch is softened, the other branches are 
almost unchanged. The agreement with the experiment 
improves compared to the "original" calculation, but the 
branch shape is clearly different, with an anti-crossing in 
the calculation that is not present in the measurements. 

It is possible to modify the nonmagnetic tetragonal 
( "original" ) calculation to agree better with the observed 




2. 'in-plane soft' 



(b) force constant (c) magnetic correlation 
'clipped' model 





FIG. 10: (Color online) Schematic of model modifications, (a) 
The "soft Fe-As" model softens both the in-plane and the out- 
of-plane components, while the "in-plane soft" reduces only 
the in-plane correlation, (b) Strength of the force constants 
between Fe atoms on an Fe-As layer in the "clipped" model. 
A strong correlation exists along the ferromagnetic direction 
as shown by the thick lines, (c) Effective magnetic exchange 
coupling constants between Fe atoms in an Fe-As layer sug- 
gested by spin wave measurement in SrFe2As2 (Ref. [6a ) and 
CaFe? As? (Re f. l66h . Strong interaction (large J) is along the 
direction of antiferromagnetic order. 



dispersion by softening only the in-plane components of 
the nearest neighbor Fe-As force constant matrix. The 
"original" model failed to fit the data, because (1) the 
in-plane polarized mode at T point (A and B in Fig. [7J 
with atomic motions A and B in Fig. [5J has an energy 
that is too high compared to the data, and (2) as one 
increases q along the <100> direction, the model pre- 
dicts an anti-crossing with a c-axis polarized mode (C 
in Fig. [7J and the atomic motion C in Fig. [5]) that is 
not observed. However, selectively reducing the in-plane 
components of the Fe-As force constant matrix by 20 % 
reduces the energy of the in-planc modes to be equal or 
just below that of the c-axis polarized mode, and avoids 
the anti-crossing, in good agreement with the measured 
data. This "in-plane soft" model [Fig. HUT a)^] is com- 
pared with the experimental data in Fig. Illf b). and, for 
being a relatively simple modification, agrees well with 
the measured dispersion. The improvement is also clearly 
seen along (3 + q 3 + q 0) in Figs.[T27a) and (b). In partic- 
ular this direction shows the necessity of softening only 
the in-plane components instead of the full force constant 
matrices as we did in Ref. H3- 

Our third modified model begins with the force con- 
stant matrices from the "magnetic" first principles cal- 
culation. To reduce the splitting between ferromagnetic 
modes (whose calculated frequencies agree with the data) 
and the antiferromagnetic ones (which are calculated 
to have a higher energy than observed), we reduce the 
<110> components of force constant matrices between 
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FIG. 11: (Color online) (a) Comparison of the experimental 
data of PrFeAsOi_y (T c = 45 K; doped-2) at Q = (3 + q 0) 
at room temperature with the "soft Fe-As" model, (b) the 
comparison with "in-plane soft" model, and (c) the compari- 
son with "clipped" model. The area of the experimental data 
as well as the line width of the calculations show the peak 
intensity on IXS spectra. See text for details. 

(a) soft Fe-As (b) in-plane soft (c) clipped 

r Mr Mr m 
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q in (3+q 3+q 0) q in (3+q ±(3+q) 0) 

FIG. 12: (Color online) As Fig. [11] but the data are from 
parent PrFeAsO and along Q = (3+q 3+q 0). In (c) "clipped" 
model the <110> and <ll0> directions are inequivalent and 
are shown in different colors (gray scales). 

nearest neighbor Fe atoms. Since the force constant be- 
tween nearest neighbor Fe atoms is much smaller (less 
than 15%) than that between nearest neighbor Fe and As, 
we "clip" these bonds by setting the force constant ma- 
trix to zero as shown in Fig. fTOT b). This "clipped" model 
is compared with the experimental data in Fig. lllf c). and 
the agreement is also relatively good. 

The "in-plane soft" and "clipped" models have reason- 
able, though not perfect, agreement with the experimen- 
tal data in other symmetry directions. For the dispersion 
along the (3 + q 0) (Fig.[H]) and (0 q 9) (not shown) the 
"clipped" model is perhaps slightly better, while for the 




Energy [meV] 

FIG. 13: (Color online) IXS spectra at Q = (3.22 3.22 -0.11) 
in parent PrFeAsO sample at 10 K around the energy trans- 
fer of 30 meV; magnification of Fig. [f^b). The symbols are 
also the same as those in Fig. [6] While the red open circles 
show the measurement at a mostly single domain region, blue 
open squares show the measurement at a multi domain region. 
Two modified models as well as the "original" first principles 
calculation are also plotted with lines. 



(3 q 0) direction (not shown) the "in-plane soft" model 
is better. However, if we consider the splitting in the 
(3 + q 3 + q 0) direction of Fig. [T^J we find that the 
"clipped" model continues to predict some splitting that 
is not observed experimentally. This is also seen in the 
comparison of the IXS spectra of the parent sample below 
Tjv at single and multi domain regions with the "original" 
ab initio calculation and two modified models in Fig. 1131 



V. DISCUSSION 

The discussion above, after extensive investigation of 
several samples of PrFeAsO \- y (y ~0, 0.1, and 0.3) and 
comparison with ab initio calculations and the modified 
models, allows the following conclusions: 

1 . There are small changes in the observed phonon in- 
tensities and dispersion when temperature or dop- 
ing is modified. 

2. The observed changes are generally much smaller 
than the rather large differences between the ob- 
served data and the ab initio calculations. 

3. Magnetic (spin-polarized LSDA) calculations tend 
to give better agreement with the measured data, 
because ferromagnetically polarized modes are soft- 
ened. However, antiferromagnetically polarized 
modes are calculated to have energies that remain 
high and do not agree with measurement. 
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4. Better agreement with the data can be obtained 
by modifying the calculation either by clipping the 
bond between antifcrromagnctically polarized iron 
atoms in spin-polarized calculations, or by soften- 
ing only the in-plane components of the force con- 
stant matrix of the Fe-As bond in nonmagnetic cal- 
culation. 

We now discuss these results in the context of other work. 
At the start we note that variations of [TH3] have already 
been reported for other iron-arsenide compounds. Our 
work then demonstrates the general applicability of these 
conclusions across many families of iron-arsenide com- 
pounds. We also have put some results on rather firmer 
footing (e.g. the lack of splitting of modes in explicitly 
twinned material). Moreover, after a precise comparison 
with "magnetic" first principles calculations, which had 
been previously considered to show good agreement with 
experimental data by others, we still find significant dif- 
ferences. Consideration of these differences leads to two 
modified models summarized in the above point |4j 

The presence of only small changes in phonon spectra 
resulting from either temperature or doping, is in agree- 
ment with work by many other authors. Phonon DOS 
measurements by either INS or IXS have been carried 
out in various iron-arsenide compounds; "1111" materi- 
als [LaFcAs(0,F) (Refs. MWR), PrFeAsOi_ y (Ref. M), 
NdFeAs(0,F) (Ref. H), Ca(Fe,Co)AsF (Ref. IH), and 
SrFeAsF (Ref. HH)], "122" materials [(Ba,K)Fe 2 As 2 
(Ref. J3l, (Ca,Na)Fe 2 As 2 (RefJS), Ca(Fe,Co) 2 As 2 
(Ref. EH), and SrFe 2 As 2 (Ref. EM and "11" system 
[FeSei-z (Ref. 0)]. Fe partial phonon DOS measure- 
ments by nuclear resonant inelastic scattering of syn - 
chrotron radiation is adapted to LaFeAs(0,F) (Ref. |67| ). 
(Ba,K)Fe 2 As 2 (Ref. HI), Ba(Fe,Co)As 2 (Ref. (H), and 
Fci+aSe (Ref. [70). Using single crystals precise phonon 
dispersions have been investigated by either IXS or INS 
in PrFcAsOi^y (Refjp), (Ba,K)Fe 2 As 2 (Ref.©, and 
SmFeAs(0,F) (Ref. l49t). All of them found only rather 
small changes, if any, with doping and temperature. The 
exception is the change in (Ba,K)Fe 2 As 2 (Ref.l44T). but it 
is observed in only one particular phonon mode. Though 
(Ca,Na)Fc 2 As 2 also shows some changes in phonon struc- 
ture with pressure^ 2 -, this system is a relatively soft mate- 
rial with a structural transition to "collapsed tetragonal" 
phase. Raman^ir— and infrared (IR) spectroscopy 8 ^—, 
or femtosecond-resolved pump-probe reflectivity 8 ^— , are 
also used to investigate phonons in these materials, and 
no significant change is observed. Our result of a gen- 
erally weak dependence of the phonon spectra on dop- 
ing and temperature is then reasonably consistent with 
previous work. Higher resolution work does show some 
changes typically at the level of 0.5 meV. 

Compared to the small temperature/doping depen- 
dences, the discrepancies between the measured data 
and the ab initio calculations are relatively large (i.e., 
a softening of observed Fe modes). Moreover, the dis- 
crepancies become smaller if magnetism is included in 
the calculations. In general, this is also in good agree- 



ment with previous results. Many phonon DOS as well 
as phonon dispersion measurements show the softening 
of the experimental data compared with nonmagnetic 
ab initio calculation by several meV on selected Fe and 
As phonon modes 3 ^ 38 ' 45-51 ' 92 . Among them some re- 
port the improvement possible using magnetic calcula- 
tions; phonon DOS of BaFe 2 As 2 (Refs.[92|(93|), i?eFe 2 As 2 
and i?eFeAsF (Ref. l5ll) . phonon dispersion of CaFe 2 As 2 
(Ref. E3), SmFeAs(0,F) (Ref. ©, (Ba,K)Fe 2 As 2 and 
Ba(Fe,Co) 2 As 2 (Refs. I4al48[ ). Here we note that the 
most pronounced softening is observed for high energy 
plane-polarized modes (AE ~30-35 meV) in our present 
work, while most of the previous results were focused to 
AE ^20 meV especially in the phonon dispersions. 

Given the generally good agreement of our calculations 
of the ferromagnetically polarized modes with the mea- 
surements on the magnetically ordered parent below Tjv, 
the remaining discrepancy between the calculated high 
energy for the antiferromagnetically polarized modes and 
the measurements which show a lower energy (essentially 
the same as the ferromagnetically polarized modes) is 
surprising. This is a strong result in that essentially all 
of our measurements only show the lower energy modes, 
and is highlighted by the complete lack of splitting (to 
the 0.5 meV level) for the occasion when we specifically 
examined a twinned portion of the sample. A similar 
lack of splitting was also observed by Reznik et a/4 8 - for 
BaFe 2 As 2 , though without explicit confirmation of twin- 
ning. 

We tried to understand the disagreement of data and 
calculations by ad hoc modifications to our calculation 
results, and arrived at two models, which have a com- 
pletely different approach, but both of which show rela- 
tively good agreement with the observed data. 

First we consider the "clipped" model. In this model 
we use the "magnetic" calculation, and include additional 
in-plane anisotropy to cancel the anisotropy originating 
from the magnetic order. The observed in-plane phonon 
dispersion of even antiferromagnetically ordered phase is 
surprisingly isotropic in the Fe-As plane. Therefore, ad- 
ditional anisotropy in the force constant matrices is indis- 
pensable, if we based it on "magnetic" calculation which 
it is now widely believed to be the reasonable calculation 
of the iron-arsenide system. The anisotropy in the Fe- 
As plane, included in "clipped" model, is also suggested 
by other experiments. For example, from magnetic exci- 
tation measurements by INS large anisotropic exchange 
coupling is suggested in CaFc 2 As 2 [SJ a — 49.9 ± 9.9 
and SJ b = -5.7 ± 4.5 (Ref. Ml or 24 < J a < 37 and 
7 < J b < 20 (Ref. H)] or in FeSeo. 5 Te .5 (Ref. HI). 
For these the strong correlation is along the ferromag- 
netic direction in the "clipped" model [Fig. nUTb 1 )] , while 
the spin-spin coupling is stronger in the antiferromag- 
netic direction [Fig. [TUtc)]. Moreover, recently orbital 
ordering in iron-pnictides has been discussed. Using 
the polarized laser angle-resolved photoemission spec- 
troscopy (ARPES) two- fold symmetry is observed in elec- 
tronic structure of the Fe id xz orbital in BaFe 2 As 2 be- 
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low Tjy (Ref. [9^). Such orbital ordering also causes in- 
plane anisotropy and even might explain 9 - 7 - the above- 
mentioned in-plane magnetic anisotropy of the spin 
waves. These anisotropic effects seem to suggest some 
support for the "clipped" model. However, the residual 
anisotropy of mode energies in the "clipped" model is 
larger than observed, and attempts to reduce it to the 
level measured (e.g., by including an attractive term) are 
not successful. Given the evidences for anisotropy (e.g., 
in spin wav o 66 ' 94 , transport measurements 9 ^, small (~0.5 
meV) splitting of phonon modes in Raman scattering 7 - 6 -) 
it seems probable that the reality is a nearly isotropic 
model, with some very small anisotropy. However, the 
reason the direct ab initio calculations so strongly over- 
estimate the anisotropy is not clear. 

The "in-plane soft" model is an alternative based on 
the "original" microscopically in-plane isotropic calcula- 
tion. We slightly soften only the in-plane components of 
force constants while preserving the out-of-plane compo- 
nent. This is an improved version of our previous "soft 
Fe-As" model^i, based on the nonmagnetic calculation. 
However, it shows as good as or even better agreement 
with the measured phonon structure compared to the 
"clipped" model. Therefore, it seems worth considering 
nonmagnetic calculations again, because experimentally 
the magnetic order does not strongly affect the phonon 
structure, and the dispersion is essentially isotropic in 
many iron-arsenides. 

Finally we consider the effect of magnetic fluctuations, 
which is widely discussed and has the possibility to re- 
duce the in-plane anisotropy. One model by Mazin and 
Johannes^ 9 -, suggested there are always antiferromagnetic 
domains, but the boundaries fluctuate. Thus depending 
on the time scale of an experimental probe relative to 
these fluctuations, different determinations may be made 
about the presence of magnetism. They suggested that 
many of the physical properties of iron-arsenide com- 
pounds, both experimental and theoretical, may be ex- 
plained by this. With this model the small observed 
Fe magnetic moment is explained as by averaging over 
zero-point and fast magnetic domain motion. Thus, 
Mossbauer experiment s 67 ' 100 and /xSRi^-iS 4 -, that probe 
relatively slow (>10 ns) timescales would see primarily 
nonmagnetic response, as they do, while probes of faster 
time scales might see magnetic response, such as the split- 
ting observed in photo-emissioni ^. INS investigations, 
in fact, show magnetic fluctuations at finite energy (^0.2 
ps time scales) even in superconducting samples and/or 
at high temperature above transitioniS 6 ."^^. 

However, it seems unlikely to us that the presence of 
fluctuations can fully explain the difference between our 
calculations and our data. Phonon time scales are given 
by the oscillation period (~0.1 ps for a 30 meV mode), 
and the phonon lifetime (~0.7 ps for a 1 meV mode) 
linewidth. Fluctuations on that are slow on the ps time 



scale should lead to splitting of modes, while those with 
timescales <0.1 ps should lead to unsplit phonons that 
see the average structure. The magnetic calculations give 
mode splittings of several meV while for the particular 
case of Fig. [3 for example, our data would support split- 
tings no larger than a few tenths of a meV. Also, the 
observed mode energy is in good agreement with the one 
of the two calculated modes, not their average. Thus, 
we speculate there is some additional neglected factor, 
not only magnetic fluctuations, that leads to the over- 
estimation of the calculated phonon anisotropy in these 
systems. 



VI. CONCLUSION 

Extensive investigation of phonon dispersion in sev- 
eral samples of PrFcAsOi_ a allow careful comparison of 
results between different samples, and with calculations. 
Small changes were visible as the material was doped into 
the superconducting state (from y = to y ~ 0.3), and 
when the parent materials was cooled below the mag- 
netic and structural phase transition. Smaller changes 
were visible when superconducting samples were cooled 
across T c . Interpretation of the observed small changes 
with doping or temperature in this system is compli- 
cated by the generally larger disagreement with calcu- 
lations. Magnetic calculations compare more favorably 
to the data, especially the fcrromagnctically polarized 
phonon modes, but show splitting between ferromagncti- 
cally and antiferromagnetically polarized modes (several 
meV) that is not at all supported by the experimental 
data (<0.5 meV). Interestingly, the in-plane isotropy of 
phonon properties can be constructed from both micro- 
scopically anisotropic (magnetic) and isotropic (nonmag- 
netic) calculations. It is clear that while the ab ini- 
tio calculations provide a reasonable rough estimate of 
phonon dispersion, they lack some fundamental ingredi- 
ent. Given the similarity of the measured response for 
both the magnetically ordered parent material below T/v 
and the other systems, however, is seems unlikely that 
the missing ingredient might be only magnetic fluctua- 
tions. 



Acknowledgments 

Work at SPring-8 was carried out under proposal 
No. 2008A2050, No. 2008B1403, No. 2009A1436, No. 
2009B1609, No. 2009B2136, and No. 2010A1296. This 
work was performed under the NIMS-RIKEN-JAEA Co- 
operative Research Program on Quantum Beam Science 
and Technology. 



12 



1 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 
J. Am. Chem. Soc. 130, 3296 (2008). 

2 C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. R. II, 
J. L. Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. L. 
Wang, et al., Nature (London) 453, 899 (2008). 

3 S. A. J. Kimber, D. N. Argyriou, F. Yokaichiya, 
K. Habicht, S. Gerischer, T. Hansen, T. Chatterji, 
R. Klingeler, C. Hess, G. Behr, et al., Phys. Rev. B 78, 
140503(R) (2008). 

4 J. Zhao, Q. Huang, C. de la Cruz, S. Li, J. W. Lynn, 
Y. Chen, M. A. Green, G. F. Chen, G. Li, Z. Li, et al., 
Nature Mater. 7, 953 (2008). 

5 Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, Z.-C. Li, 

G. -C. Che, X.-L. Dong, L.-L. Sun, F. Zhou, et al., Chin. 
Phys. Lett. 25, 2215 (2008). 

6 H. Kito, H. Eisaki, and A. Iyo, J. Phys. Soc. Jpn. 77, 
063707 (2008). 

7 G. Wu, Y. L. Xie, H. Chen, M. Zhong, R. H. Liu, B. C. 
Shi, Q. J. Li, X. F. Wang, T. Wu, Y. J. Yan, et al., J. 
Phys.: Condens. Matter 21, 142203 (2009). 

8 M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 

9 K. Sasmal, B. Lv, B. Lorenz, A. M. Guloy, F. Chen, Y.- 
Y. Xue, and C.-W. Chu, Phys. Rev. Lett. 101, 107007 
(2008). 

10 X. C. Wang, Q. Q. Liu, Y. X. Lv, W. B. Gao, L. X. 
Yang, R. C. Yu, F. Y. Li, and C. Q. Jin, Solid State 
Communications 148, 538 (2008). 

11 D. R. Parker, M. J. Pitcher, P. J. Baker, I. Franke, T. Lan- 
caster, S. J. Blundell, and S. J. Clarke, Chemical Com- 
munications 2009, 2189 (2009). 

12 F.-C. Hsu, J.-Y. Luo, K.-W. Yeh, T.-K. Chen, T.-W. 
Huang, P. M. Wu, Y.-C. Lee, Y.-L. Huang, Y.-Y. Chu, 
D.-C. Yan, et al., Proceedings of the National Academy 
of Sciences of the United States of America 105, 14262 
(2008). 

13 Y. Mizuguchi, F. Tomioka, S. Tsuda, T. Yamaguchi, and 
Y. Takano, Appl. Phys. Lett. 93, 152505 (2008). 

14 H. Ogino, Y. Matsumura, Y. Katsura, K. Ushiyama, 
S. Horii, K. Kishio, and J. Shimoyama, Superconductor 
Science and Technology 22, 075008 (2009). 

15 X. Zhu, F. Han, G. Mu, P. Cheng, B. Shen, B. Zeng, and 

H. -H. Wen, Phys. Rev. B 79, 220512(R) (2009). 

16 N. Kawaguchi, H. Ogino, Y. Shimizu, K. Kishio, and 
J. Shimoyama, Applied Physics Express 3, 063102 (2010). 
H. Ogino, S. Sato, K. Kishio, J. Shimoyama, T. Tohei, 
and Y. Ikuhara, Appl. Phys. Lett. 97, 072506 (2010). 

18 H. Ogino, K. Machida, A. Yamamoto, K. Kishio, J. Shi- 
moyama, T. Tohei, and Y. Ikuhara, Superconductor Sci- 
ence and Technology 23, 115005 (2010). 

19 P. M. Shirage, K. Kihou, C.-H. Lee, H. Kito, H. Eisaki, 
and A. Iyo, Appl. Phys. Lett. 97, 172506 (2010). 

20 C.-H. Lee, A. Iyo, H. Eisaki, H. Kito, M. T. Fernandez- 
Diaz, T. Ito, K. Kihou, H. Matsuhata, M. Braden, and 
K. Yamada, J. Phys. Soc. Jpn. 77, 083704 (2008). 

21 M. J. Han, Q. Yin, W. E. Pickett, and S. Y. Savrasov, 
Phys. Rev. Lett. 102, 107003 (2009). 

22 K. Kuroki, H. Usui, S. Onari, R. Arita, and H. Aoki, Phys. 
Rev. B 79, 224511 (2009). 

23 K. Horigane, H. Hiraka, and K. Ohyama, J. Phys. Soc. 
Jpn. 78, 074718 (2009). 



24 C. de la Cruz, W. Z. Hu, S. Li, Q. Huang, J. W. Lynn, 
M. A. Green, G. F. Chen, N. L. Wang, H. A. Mook, Q. Si, 
et al., Phys. Rev. Lett. 104, 017204 (2010). 

25 T. Nomura, Y. Inoue, S. Matsuishi, M. Hirano, J. E. Kim, 
K. Kato, M. Takata, and H. Hosono, Superconductor Sci- 
ence and Technology 22, 055016 (2009). 

26 H. Takahashi, K. Igawa, K. Arii, Y. Kamihara, M. Hirano, 
and H. Hosono, Nature (London) 453, 376 (2008). 

27 S. Medvedev, T. M. McQueen, I. A. Troyan, T. Palasyuk, 
M. I. Eremets, R. J. Cava, S. Naghavi, F. Casper, 
V. Ksenofontov, G. Wortmann, et al., Nature Mater. 8, 
630 (2009). 

28 D. J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 
(2008). 

29 L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. Rev. 
Lett. 101, 026403 (2008). 

30 K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 

100, 226402 (2008). 

31 I. I. Mazin, M. D. Johannes, L. Boeri, K. Koepernik, and 
D. J. Singh, Phys. Rev. B 78, 085104 (2008). 

32 T. Yildirim, Phys. Rev. Lett. 102, 037003 (2009). 

33 L. Boeri, M. Calandra, I. I. Mazin, O. V. Dolgov, and 
F. Mauri, Phys. Rev. B 82, 020506(R) (2010). 

34 F. Yndurain, Europhys. Lett. 94, 37001 (2011). 

35 Y. Qiu, M. Kofu, W. Bao, S.-H. Lee, Q. Huang, 
T. Yildirim, J. R. D. Copley, J. W. Lynn, T. Wu, G. Wu, 
et al., Phys. Rev. B 78, 052508 (2008). 

36 A. D. Christianson, M. D. Lumsden, O. Delaire, M. B. 
Stone, D. L. Abernathy, M. A. McGuire, A. S. Sefat, 
R. Jin, B. C. Sales, D. Mandrus, et al., Phys. Rev. Lett. 

101, 157004 (2008). 

37 T. Fukuda, A. Q. R. Baron, S. Shamoto, M. Ishikado, 
H. Nakamura, M. Machida, H. Uchiyama, S. Tsutsui, 
A. Iyo, H. Kito, et al., J. Phys. Soc. Jpn. 77, 103715 
(2008). 

38 M. L. Tacon, M. Krisch, A. Bosak, J.-W. G. Bos, and 
S. Margadonna, Phys. Rev. B 78, 140505(R) (2008). 

39 R. Mittal, Y. Su, S. Rols, T. Chatterji, S. L. Chaplot, 
H. Schober, M. Rotter, D. Johrendt, and T. Brueckel, 
Phys. Rev. B 78, 104514 (2008). 

40 R. Mittal, Y. Su, S. Rols, M. Tegel, S. L. Chaplot, 
H. Schober, T. Chatterji, D. Johrendt, and T. Brueckel, 
Phys. Rev. B 78, 224518 (2008). 

41 R. Mittal, S. Rols, M. Zbiri, Y. Su, H. Schober, S. L. 
Chaplot, M. Johnson, M. Tegel, T. Chatterji, S. Matsu- 
ishi, et al., Phys. Rev. B 79, 144516 (2009). 

42 R. Mittal, L. Pintschovius, D. Lamago, R. Heid, K.-P. 
Bohnen, D. Reznik, S. L. Chaplot, Y. Su, N. Kumar, S. K. 
Dhar, et al., Phys. Rev. Lett. 102, 217001 (2009). 

43 D. Phelan, J. N. Millican, E. L. Thomas, J. B. Lefio, 
Y. Qiu, and R. Paul, Phys. Rev. B 79, 014519 (2009). 

44 C.-H. Lee, K. Kihou, K. Horigane, S. Tsutsui, T. Fukuda, 
H. Eisaki, A. Iyo, H. Yamaguchi, A. Q. R. Baron, 
M. Braden, et al., J. Phys. Soc. Jpn. 79, 014714 (2010). 

45 D. Reznik, K. Lokshin, D. C. Mitchell, D. Parshall, 
W. Dmowski, D. Lamago, R. Heid, K.-P. Bohnen, A. S. 
Sefat, M. A. McGuire, et al., arXiv: 0810.4941. 

46 R. Mittal, M. Zbiri, S. Rols, Y. Su, Y. Xiao, H. Schober, 
S. L. Chaplot, M. Johnson, T. Chatterji, S. Matsuishi, 
et al., Phys. Rev. B 79, 214514 (2009). 

47 S. E. Hahn, Y. Lee, N. Ni, P. C. Canfield, A. I. Goldman, 



13 



R. J. McQueeney, B. N. Harmon, A. Alatas, B. M. Leu, 
E. E. Alp, et al., Phys. Rev. B 79, 220511(R) (2009). 
D. Reznik, K. Lokshin, D. C. Mitchell, D. Parshall, 
W. Dmowski, D. Lamago, R. Heid, K.-P. Bohnen, A. S. 
Sefat, M. A. McGuire, et al., Phys. Rev. B 80, 214534 
(2009). 

M. L. Tacon, T. R. Forrest, C. Riiegg, A. Bosak, A. C. 
Walters, R. Mittal, H. M. R0nnow, N. D. Zhigadlo, S. Ka- 
trych, J. Karpinski, et al., Phys. Rev. B 80, 220504(R) 

(2009) . 

M. L. Tacon, T. R. Forrest, C. Riiegg, A. Bosak, J. Noff- 

singer, A. C. Walters, P. Toulemonde, A. Palenzona, N. D. 

Zhigadlo, J. Karpinski, et al., arXiv: 1005.4165. 

M. Zbiri, R. Mittal, S. Rols, Y. Su, Y. Xiao, H. Schober, 

S. L. Chaplot, M. R. Johnson, T. Chatterji, Y. Inoue, 

et al., J. Phys.: Condens. Matter 22, 315701 (2010). 

M. Ishikado, S. Shamoto, H. Kito, A. Iyo, H. Eisaki, 

T. Ito, and Y. Tomioka, Physica C 469, 901 (2009). 

C. Hess, A. Kondrat, A. Narduzzo, J. E. Hamann-Borrero, 

R. Klingeler, J. Werner, G. Behr, and B. Biichner, Euro- 

phys. Lett. 87, 17005 (2009). 

M. Rotter, M. Tegel, D. Johrendt, I. Schellenberg, 
W. Hermes, and R. Pottgen, Phys. Rev. B 78, 020503(R) 
(2008). 

J. Zhao, Q. Huang, C. de la Cruz, J. W. Lynn, M. D. 
Lumsden, Z. A. Ren, J. Yang, X. Shen, X. Dong, Z. Zhao, 
et al., Phys. Rev. B 78, 132504 (2008). 
A. Q. R. Baron, Y. Tanaka, S. Goto, K. Takeshita, 
T. Matsushita, and T. Ishikawa, J. Phys. Chem. Solids 
61, 461 (2000). 

A. Q. R. Baron, J. P. Sutter, S. Tsutsui, H. Uchiyama, 
T. Masui, S. Tajima, R. Heid, and K.-P. Bohnen, J. Phys. 
Chem. Solids 69, 3100 (2008). 

G. Kresse and J. Hafner, Phys. Rev. B 47, (R)558 (1993). 
G. Kresse and J. Furthmiiller, Computational Materials 
Science 6, 15 (1996). 

G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 
78, 4063 (1997). 

S. Ishibashi, K. Terakura, and H. Hosono, J. Phys. Soc. 
Jpn. 77, 053709 (2008). 

J. Zhao, D.-X. Yao, S. Li, T. Hong, Y. Chen, S. Chang, 
W. R. II, J. W. Lynn, H. A. Mook, G. F. Chen, et al., 
Phys. Rev. Lett. 101, 167203 (2008). 
J. Zhao, D. T. Adroja, D.-X. Yao, R. Bewley, S. Li, X. F. 
Wang, G. Wu, X. H. Chen, J. Hu, and P. Dai, Nature 
Physics 5, 555 (2009). 

S. Higashitaniguchi, M. Seto, S. Kitao, Y. Kobayashi, 
M. Saito, R. Masuda, T. Mitsui, Y. Yoda, Y. Kamihara, 
M. Hirano, et al., Phys. Rev. B 78, 174507 (2008). 
S. Tsutsui, C.-H. Lee, C. Tassel, Y. Yoshida, Y. Yoda, 
K. Kihou, A. Iyo, and H. Eisaki, J. Phys. Soc. Jpn. 79, 
013706 (2010). 

O. Delaire, M. S. Lucas, A. M. dos Santos, A. Subedi, 
A. S. Sefat, M. A. McGuire, L. Mauger, J. A. Munoz, 
C. A. Tulk, Y. Xiao, et al., Phys. Rev. B 81, 094504 

(2010) . 

V. Ksenofontov, G. Wortmann, A. I. Chumakov, T. Gasi, 
S. Medvedev, T. M. McQueen, R. J. Cava, and C. Felser, 
Phys. Rev. B 81, 184510 (2010). 

A. P. Litvinchuk, V. G. Hadjiev, M. N. Iliev, B. Lv, 



A. M. Guloy, and C. W. Chu, Phys. Rev. B 78, 060503(R) 
(2008). 

S. C. Zhao, D. Hou, Y. Wu, T. L. Xia, A. M. Zhang, G. F. 

Chen, J. L. Luo, N. L. Wang, J. H. Wei, Z. Y. Lu, et al., 

Supercond. Sci. Technol. 22, 015017 (2009). 

Y. Gallais, A. Sacuto, M. Cazayous, P. Cheng, L. Fang, 

and H. H. Wen, Phys. Rev. B 78, 132509 (2008). 

L. Zhang, T. Fujita, F. Chen, D. L. Feng, S. Maekawa, 

and M. W. Chen, Phys. Rev. B 79, 052507 (2009). 

M. Rahlenbeck, G. L. Sun, D. L. Sun, C. T. Lin, 

B. Keimer, and C. Ulrich, Phys. Rev. B 80, 064509 (2009). 
L. Chauviere, Y. Gallais, M. Cazayous, A. Sacuto, M. A. 
Measson, D. Colson, and A. Forget, Phys. Rev. B 80, 
094504 (2009). 

K.-Y. Choi, D. Wulferding, P. Lemmens, N. Ni, S. L. 
Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 212503 
(2008). 

K.-Y. Choi, P. Lemmens, I. Eremin, G. Zwicknagl, 
H. Berger, G. L. Sun, D. L. Sun, and C. T. Lin, J. Phys.: 
Condens. Matter 22, 115802 (2010). 

M. Granath, J. Bielecki, J. Holmlund, and L. Borjesson, 
Phys. Rev. B 79, 235103 (2009). 

P. Kumar, A. Kumar, S. Saha, D. V. S. Muthu, 
J. Prakash, S. Patnaik, U. V. Waghmare, A. K. Ganguli, 
and A. K. Sood, Solid State Communications 150, 557 
(2010). 

T.-L. Xia, D. Hou, S. C. Zhao, A. M. Zhang, G. F. Chen, 
J. L. Luo, N. L. Wang, J. H. Wei, Z.-Y. Lu, and Q. M. 
Zhang, Phys. Rev. B 79, 140510(R) (2009). 
K. Okazaki, S. Sugai, S. Niitaka, and H. Takagi, Phys. 
Rev. B 83, 035103 (2011). 

A. M. Zhang, J. H. Xiao, J. B. He, D. M. Wang, G. F. 
Chen, and Q. M. Zhang, arXiv: 1101.2168. 
S. Sugai, Y. Mizuno, K. Kiho, M. Nakajima, C. H. Lee, 
A. Iyo, H. Eisaki, and S. Uchida, Phys. Rev. B 82, 
140504(R) (2010). 

S. I. Mirzaei, V. Guritanu, A. B. Kuzmenko, C. Senatore, 
D. van der Marel, G. Wu, R. H. Liu, and X. H. Chen, 
arXiv: 0806.2303. 

A. Akrap, J. J. Tu, L. J. Li, G. H. Cao, Z. A. Xu, and 

C. C. Homes, Phys. Rev. B 80, 180502(R) (2009). 

D. Wu, N. Barisic, N. Drichko, S. Kaiser, A. Faridian, 
M. Dressel, S. Jiang, Z. Ren, L. J. Li, G. H. Cao, et al., 
Phys. Rev. B 79, 155103 (2009). 

E. van Heumen, Y. Huang, S. de Jong, A. B. Kuzmenko, 
M. S. Golden, and D. van der Marel, Europhys. Lett. 90, 
37005 (2010). 

B. Mansart, D. Boschetto, A. Savoia, F. Rullier-Albenque, 
A. Forget, D. Colson, A. Rousse, and M. Marsi, Phys. 
Rev. B 80, 172504 (2009). 

T. Mertelj, P. Kusar, V. V. Kabanov, L. Stojchevska, 
N. D. Zhigadlo, S. Katrych, Z. Bukowski, J. Karpinski, 
S. Weyeneth, and D. Mihailovic, Phys. Rev. B 81, 224504 
(2010). 

H. Takahashi, Y. Kamihara, H. Koguchi, T. Atou, 
H. Hosono, I. Katayama, J. Takeda, M. Kitajima, and 
K. G. Nakamura, J. Phys. Soc. Jpn. 80, 013707 (2011). 
M. Zbiri, H. Schober, M. R. Johnson, S. Rols, R. Mittal, 
Y. Su, M. Rotter, and D. Johrendt, Phys. Rev. B 79, 
064511 (2009). 

T. Yildirim, Physica C 469, 425 (2009). 
S. O. Diallo, V. P. Antropov, T. G. Perring, C. Bro- 
holm, J. J. Pulikkotil, N. Ni, S. L. Bud'ko, P. C. Canfield, 
A. Kreyssig, A. I. Goldman, et al., Phys. Rev. Lett. 102, 



14 



187206 (2009). 

5 S.-H. Lee, G. Xu, W. Ku, J. S. Wen, C. C. Lee, 
N. Katayama, Z. J. Xu, S. Ji, Z. W. Lin, G. D. Gu, et al., 
Phys. Rev. B 81, 220502(R) (2010). 

6 T. Shimojima, K. Ishizaka, Y. Ishida, N. Katayama, 
K. Ohgushi, T. Kiss, M. Okawa, T. Togashi, X.-Y. Wang, 
C.-T. Chen, et al., Phys. Rev. Lett. 104, 057002 (2010). 

7 W. Lv, F. Kriiger, and P. Phillips, Phys. Rev. B 82, 
045125 (2010). 

8 J.-H. Chu, J. G. Analytis, K. D. Greve, P. L. McMahon, 
Z. Islam, Y. Yamamoto, and I. R. Fisher, Science 329, 
824 (2010). 

9 I. I. Mazin and M. D. Johannes, Nature Physics 5, 141 
(2009). 

S. Kitao, Y. Kobayashi, S. Higashitaniguchi, M. Saito, 
Y. Kamihara, M. Hirano, T. Mitsui, H. Hosono, and 
M. Seto, J. Phys. Soc. Jpn. 77, 103706 (2008). 

1 H. Luetkens, H.-H. Klauss, R. Khasanov, A. Amato, 
R. Klingeler, I. Hellmann, N. Leps, A. Kondrat, C. Hess, 
A. Kohler, et al., Phys. Rev. Lett. 101, 097009 (2008). 

2 J. P. Carlo, Y. J. Uemura, T. Goko, G. J. MacDougall, 
J. A. Rodriguez, W. Yu, G. M. Luke, P. Dai, N. Shannon, 



S. Miyasaka, et al., Phys. Rev. Lett. 102, 087001 (2009). 
S. Takeshita, R. Kadono, M. Hiraishi, M. Miyazaki, 
A. Koda, S. Matsuishi, and H. Hosono, Phys. Rev. Lett. 
103, 027002 (2009). 

K. Ohishi, Y. Ishii, K. Miyazawa, H. Fukazawa, I. Watan- 

abe, Y. Kohori, P. M. Shirage, H. Kito, A. Iyo, and 

H. Eisaki, J. Phys. Soc. Jpn. 80, 024703 (2011). 

F. Bondino, E. Magnano, M. Malvestuto, F. Parmigiani, 

M. A. McGuire, A. S. Sefat, B. C. Sales, R. Jin, D. Man- 

drus, E. W. Plummer, et al., Phys. Rev. Lett. 101, 267001 

(2008). 

A. D. Christianson, E. A. Goremychkin, R. Osborn, 
S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. 
Todorov, H. Claus, D. Y. Chung, M. G. Kanatzidis, et al., 
Nature (London) 456, 930 (2008). 

M. Ishikado, R. Kajimoto, S. Shamoto, M. Arai, A. Iyo, 
K. Miyazawa, P. M. Shirage, H. Kito, H. Eisaki, S. Kim, 
et al., J. Phys. Soc. Jpn. 78, 043705 (2009). 
S. Wakimoto, K. Kodama, M. Ishikado, M. Matsuda, 
R. Kajimoto, M. Arai, K. Kakurai, F. Esaka, A. Iyo, 
H. Kito, et al., J. Phys. Soc. Jpn. 79, 074715 (2010). 



